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Abstract 

In low temperature limit, we study electron counting statistics of a disordered 

conductor. We derive an expression for the distribution of charge transmitted 

over a finite time interval by using a result from the random matrix theory of 

quasi one dimensional disordered conductors. In the metallic regime, we find 

that the peak of the distribution is Gaussian and shows negligible sample to 

sample variations. We also find that the tails of the distribution are neither 

Gaussian nor Poisson and exhibit strong sample to sample variations. 
PACS numbers: 72.10.Bg, 73.50.Fq, 73.50.Td 



1 



Typeset using REVTf^X 



The physics of current fluctuations at low temperature presents an interesting quantum 
mechanical problem. Classical Johnson-Nyquist noise formula [jT| gives a good description 
of current fluctuations due to thermal fluctuations. However, at low temperature thermal 
fluctuations are small and a new type of noise becomes important. At low temperature, the 
quantum nature of the current and the discreteness of electron charge is the main source of 
current fluctuations and due to these reasons, this noise is called "quantum shot noise" . 

Lots of the low temperature current fluctuation studies deal with a disordered conductor 
because it has a simple and well established mathematical description based on Landauer's 
approach ||. Lesovik H and Yurke and Kochanski || studied the quantum shot noise in a 
two-terminal conductor using this approach and found an expression for noise power which 
is a factor 1 — T off the classical shot noise, where T is a transmission coefficient. This 
analysis was generalized to a multiterminal conductor by Biittiker ||, and he also found the 
reduction of the noise. Physically, the noise reduction is due to the Fermi statistics that 
leads to correlation of transmission events. 

In Landauer's approach, details of current transport are determined by transmission co- 
efficients and there have been many works on the distribution of the coefficients. In the past 
decade, the random matrix theory of disordered conductors was developed 0, motivated 
by theoretical discovery and experimental observation of the universal conductance fluctu- 
ations. The theory succeeded in providing a complete characterization of the distribution 
and it also succeeded in providing the insight on the origin of the universal conductance 
fluctuations. One of the fundamental results in the random matrix theory is the universal- 
ity of the distribution in the metallic regime and the universality provides a link between 
the microscopically calculated transmission coefficients and the macroscopically measurable 
conductance. 

Application of the universal distribution to current fluctuations also provides insights on 
the current fluctuations of a disordered conductor. Beenakker and Biittiker |7j calculated 
the sample averaged noise power using the universal distribution and found that it depends 
only on the conductance and that it is one-third of the classical value. 
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Noise power is a good measure of noise magnitude and its study revealed the reduction 
of noise due to Fermi correlation. However, compared to Johnson-Nyquist noise, our un- 
derstanding of quantum shot noise is limited and not many things are known besides the 
noise magnitude. Our goal in this paper is to explore the physics of low temperature cur- 
rent fluctuations beyond the scope of noise power. For our purpose, it is useful to look at 
the behavior of current fluctuations in the time domain, which brings one to the notion of 
counting statistics of charge transmitted in a conductor over fixed time. A previous study of 
the counting statistics for a single channel conductor revealed that the attempts to transmit 
electrons are highly correlated and almost periodic in time, which leads to binomial statistics 



In the time domain study of current fluctuations, the charge Q(t) (measured in the units 
of e) transmitted over a time interval t is the quantity of interest and the probability dis- 
tribution P(Q{t)) reflects the nature of the current fluctuations. Even at zero temperature, 
P{Q{t)) has finite peak width due to the quantum nature of current. To study P(Q(t)), it 
is useful to introduce the characteristic function 



because in many cases, x(A) is easier to calculate than P(Q(t)) itself. x(A) is a Fourier 
transform of P(Q(t)) and so once xW is obtained, we can either take an inverse Fourier 
transform of it to get an explicit expression for P(Q(t)), or expand it as follows to get 
cumulants of the distribution: 



fc=i "'■ 

In the linear transport regime, we derive a general expression for xW i n terms of trans- 
mission coefficients and by combining it with the transmission coefficients distribution for 
quasi one dimensional conductors, we show 



X(\) = e iQX P(Q) for - tt < A < tt , 



(1) 



integer Q 




(2) 




(3) 
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where V is the dc voltage, G = g(e 2 /h) = (Nl/L)(e 2 /h) is the average conductance, and 
the bar on the left hand side represents the sample average. Cumulant expansion of Eq. (|3|) 
implies that in average, for GVt/e ^> 1, P(Q(t)) has a Gaussian peak at GVt/e with 
((Q 2 {t))) = GVt/3e. It also implies that even though the peak is Gaussian, the tails show 
deviation from both Gaussian and Poisson distributions. We estimate sample to sample 
variations of P(Q(t)) by studying variances of various quantities and find that for GVt/e ^> 1 
sample to sample variations of P(Q(t)) appear only in the tails of P(Q(t)) and that around 
the peak, P(Q{t)) is universal. 

Before we present the derivation of above result, we stress that there are two kinds of 
averages involved. To avoid confusion, we will use a bar( TTT ) for an ensemble average, or an 
average over samples, and a bracket ((• • •)) for a quantum average, or a quantum expectation 
value. Also we reserve a double bracket (((• • •))) for a cumulant of a quantum expectation 
value and "var"(var(- • •)) for • • - 2 — ttt2 . 

Now, let us derive Eq. (|3|). Following the Landauer's approach 0, we consider a con- 
ductor sandwiched between two perfect leads. In the linear transport regime, scattering 
properties of a conductor are described by a unitary scattering matrix S that relates incom- 
ing and outgoing amplitudes, Il(r) and Ol(r)'- 



<o L \ 



(4) 



S 

where the subscripts L and R stand for the left and the right leads. 

The unitarity of S is due to the current conservation, and it allows a system to be 
decomposed into independent channels 0. Then the decomposition motivates one to study 
single channel transport first, where a transmission coefficient T determines the transport. 
Recently, the counting statistics of the single channel transport was studied || . In the low 
temperature limit (&bT eV), the characteristic function Xi(^) °f a single channel system 
becomes 

Xi(A) = (pe* x + q) M , (5) 



where p = T, q = 1 — T, M = eVt/h [JTO]] and M > 1 is assumed. The inverse Fourier 
transform of Eq. (|5]) gives the binomial distribution, which implies that the intervals between 
subsequent attempts to transmit electrons are quite regular. This regularity is due to Pauli 
exclusion principle. 

Having the characteristic function of a single channel, we write the total characteristic 
function %(A) product, 

X(A) = nC^e* + 1 - 2})", (6) 
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where Tj is a transmission coefficient of channel j. The product form Eq. (pj) follows from 
the mutual independence of channels. By taking logarithm of Eq. (pj), we get 



log X (A) = M J2 log (t^ + 1 - T^j . 
and by expanding Eq. (|7|) in terms of A, we find 



(7) 
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((^)))=ME(T(1-T)- T\ t=Tj . (8) 
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We note that both logx(A) and {{Q (t)}} are linear statistics of T^'s. 

Current fluctuations are determined by distribution of transmission coefficients and the 
distribution varies from sample to sample even though samples have the same macroscopic 
parameters. Therefore in principle each sample exhibits distinctive current fluctuations. 
However according to the random matrix theory of disordered conductors, in the metallic 
regime(l <C g <C N) where N is the number of channels, the distribution approaches a 
universal one ||. This result provides a motivation to approximate the sample-dependent 
distribution by the universal one. To exploit the universal distribution, we introduce new 
variables Uj's and the density function D{y) defined by Tj = 1/ cosh 2 uj and D(u)du = 
D(T)dT, where D{T) is the density function of 2}'s. According to Ref. 0, D(v) is uniform 
over a wide range of u, 

D{y) = g for v < u c . (9) 
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We combine Eq. (0) with the universal distribution to obtain 

log X (A)=Q 0j / o ^logl—^ + lj, (10) 



where Qq = gM. In Eq. (10) the upper limit v c is replaced by infinity, which is valid in the 
metallic regime because for large the integrand is exponentially small. The evaluation of 
the integral then leads to Eq. (^|). We note that because logx(A) is a linear statistic, the 
universal distribution approximation is equivalent to taking an average over samples. 

Cumulants are useful in understanding features of the probability distribution. By using 
the formula Eq. (0), we obtain sample averaged cumulants 



mm 


— Qo, 




— 


((Q 3 (t))) 


= j^Qo, 




= _ 105^°' 








= 231^0, 




- 27 n 

5005^°' 


((Q 8 (t)}} 


= -7T5Q0, 




233 /n 
36465^°' 


((Q w (t))) 


6823 q 
969969^° 



(11) 



The first cumulant is trivial. It is just a definition of G and it shows where the peak of 
P(Q(t)) is. The second cumulant measures the width square of the peak. It is also directly 
related to the noise power P = Jdt ((/(0)/(t))), a widely used measure of noise magnitude, 
by ((Q 2 (t))) = tP for large t, and its ensemble average is one-third of the classical value 
Qo, as first pointed out by Beenakker and Buttiker 0. The third and the fourth cumulants 
are measures of skewness and sharpness of the peak, respectively and they are related to 
3 and 4-point current current correlation functions by similar relations. We note that all 



cumulants are proportional to Qo and that for Q 1, ((Q(t))) ^> ((Q k {t)))- Therefore 
the peak of the distribution P(Q(t)) is Gaussian for large conductance limit or long time 
limit. This result is quite expected from the central limit theorem. Now to see the tails of 
P(Q(t)), we study higher order cumulants. From Eq. (|10D, we obtain a general formula for 
the ensemble averaged k-th order cumulants 



Qo r , 1 r j_ ^ q 



fc-i 



((Q k (t)}} = -i k ^- / dx , \ dqe~ iqx ^ -, (12) 

U ^ W// 4 i-» Vl + e-^i-oo y sinh(vrg - z0+) V ; 
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and by using the steepest descent method twice, we obtain the asymptotics for large k, 
Q n (k — lV (—1)^ for even k , 

((Q k m - J$=i { ^ k+1 (13) 

( 2?r ) Vfc I ( _!)^ for odd fc . 

The high order cumulants diverge as factorial, which suggests that at the tails, P(Q(t)) is 
different from both Gaussian distribution and Poisson distribution which describes the clas- 



sical current fluctuations. In comparison, ((Q k (t)}} = for k > 2 for Gaussian distribution 



and {(Q k (t)}) = Qo for k > 1 for Poisson distribution. 

It is known that in the presence of time reversal symmetry, there are order M corrections 



to ((Q(t))} and ((Q 2 (t))) due to weak localization [11], and it is legitimate to expect the 
same kind of corrections to higher order cumulants. However, because we are interested in 
the metallic regime, these corrections are small by a factor g and we will ignore them. 

A proper next step is to estimate the magnitude of sample to sample variations of 
P(Q(t)). Here instead of logx(A), we examine variance of ((Q k {t))) to see the variations of 
P(Q(t)). ((Q k {t))) is a linear statistic and the general formula for the variance of a linear 
statistic A = J2j a {Tj) is obtained recently by Beenakker and Rejaei [ 12] , and Chalker and 
Mecedo @, 

i as 1 f°° ii ka 2 (k) . 
var (A) = — - / dk- — , 14 

^ ' p7l 2 Jo 1 + COth(iTTfc) V ' 

a(k) — 2 / dv a ( § — ^jcos/cz/, (15) 

Jo Vcosh v) 

where (3 is a symmetry constant, 1, 2, or 4 depending on the symmetry. We use this formula 
to obtain 

var(«Q(t)») = ^M 2 , var («Q 2 (t)») = ^M 2 , var («Q 3 (t)») = Ti ^M 2 , • • • . 

(16) 

We note that for low order cumulants, ((Q k (t)}) is larger than var ( y ((Q k (t))}^j by at least 
a factor of g 2 , which is large in the metallic regime. Low order cumulants decide the shape 
of P(Q(t)) around the peak and therefore the small variance of low order cumulants implies 
that the peak shape shows little sample to sample variations, that is, it is almost universal. 



To see the behavior of higher order cumulants, we obtain an asymptotic form of the 



variance from an approximate variance formula in Ref. jn 
According to Eq. (|T7|), for high order cumulants, var (^((Q k (t)))^j becomes larger than 



{(Q k (t)}} due to its rapidly growing factorial factor, which suggests that the tails of P(Q(t)) 
show large sample to sample variations. We argue that this rapid growth of var (((Q fe (t)))) is 
not an artifact of the approximate variance formula used above because it assumes stronger 
spectral rigidity than the formula Eq. (|l4l , |l5|) and it has a tendency to slightly underesti- 
mate variances. Therefore large sample to sample variations at the tails of P(Q(t)) obtained 
above is not an artifact of the approximation. 



In the above, we derived the shape of P(Q(t)) by examining logx(A) and its cumulant 



expansion instead of logx(A), which might be an intuitively more appropriate ensemble 



average because it is directly related to P(Q(t)). However we argue that in contrast to the 



intuition, logx(A) is an appropriate ensemble average for the study of current fluctuations. 
One reason is that as we remarked earlier, a fc-point current current correlation function is 



linearly related to ({Q (£))), whose ensemble average can be obtained from logx(A) by a 



simple expansion. Another reason is that as we show later, logx(A) either becomes identical 



to logx(A) at short time limit, or is dominated by the conductance fluctuations instead of 
the current fluctuations. 



Calculation of x(A) is not simple because x(A) is not a linear statistic. Muttalib and Chen 



15] did this calculation recently by large N limit continuum approximation and showed that 



at long time limit, logx(A) becomes quite different from logx(A). Here we present improved 
calculation by a perturbation method and we believe that our calculation clarifies the reason 
why two averages become quite different at long time limit. 

Because xW is n °t a linear statistic, we need joint probability distribution of transmis- 
sion coefficients to average it over ensembles. After standard variable change, T = 1/(1 + x), 
the joint probability distribution P({x}) is 
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P({x}) = exp 1(3 Yl V(x a , x b )+(3Y, U(x a ) . (18) 

\ a<b a J 

We choose V(x,y) = (l/2)log(x — y) + (1/2) log^rcsinh 2 -^ — arcsinh 2 ^^) and U(x) = 
g arcsinh 2 (y / a;) based on the exact calculation of the joint probability distribution function 



for (3 = 2 by Beenakker and Rejaei | 12| . Then 



x(A) = |f, 

Z M =IUa dx a exp (ft Ea<b V(x a , X b ) +f3Ea U(x a ) + £ a M log , (19) 

Z = JUa dx a exp (ft Y,a<b V(x a , X h ) + PY.aU (x a )) . 



By expanding logx(A) in terms of M, we find 



logx(A) = log X (A) +var (log X (A)) + 0(M 3 ), (20) 



and from the formula Eq. ( fl4] , [l5"l) , we obtain 



M 2 / , arcsinhVe^ - 1 1 



logx(A) = ^MarcsinhVe^ - 1 - — 3 log ^== + ~ lX + 0(M 3 ). (21) 

P V V e lX - 1 2 I 



Note that for M -C g(short time limit), logx(A) reduces to logx(A). We expand logx(A) in 



terms of A to see features of P(Q(t)): 
log^A) = 9 M(«) + (f M + 1m=) <f + (l« + JLm* + £>(*•)) <f + . 



(22) 



The first expansion coefficient shows that (Q(t)) = Qo = gM, which is trivial. The second 



expansion coefficient, (Q 2 (t)} - (Q(t)) = ((Q 2 (t))) + var((Q(t)» = Q /3 + (2/15/3)M 2 



indicates that the peak width of P(Q(t)) has two contributions. The first contribution is 
related to the noise power, and the second one to the universal conductance fluctuations 
because v&r((Q(t))) is proportional to the variance of the conductance. (The factor 2/15/3 is 
precisely the variance of the dimensionless conductance.) Note that as t — > oo, the second 
contribution becomes dominant over the first one. It can be shown that the k-th order 
expansion coefficient contains k different contributions and at long time limit, the most 
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dominant contribution, which is proportional to M k , is related to the k-th cumulant of the 
conductance fluctuations. From this analysis we see that the behavior of x(A) for large t is 
governed by the conductance fluctuations instead of the current fluctuations. 

As a short remark, we report that the continuum approximation calculation of logx(A), 
as suggested by Muttalib and Chen JOJ, produces the same result as Eq. ( plj ) up to M 2 
order. 

In summary, we examine the counting statistics of charge to study the low temperature 
current fluctuations. By calculating the characteristic function of the probability distribution 
P(Q(t)), we find that P(Q(t)) has a Gaussian peak at Qo with ((Q(t))) = Qo/3 and we also 
find that the tails of P(Q(t)) are different from the tails of Gaussian and classical Poisson 
distribution. By studying the variances of the cumulants, we establish that even though 
the peak location of P(Q(t)) varies from sample to sample due to universal conductance 
fluctuations, the peak shape of P(Q(t)) is universal in the metallic regime, and that the 
sample to sample variations show up only at the tails of P(Q(t)). 
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